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Abstract 



Using the principles of the conformal quantum field theory and the finite size 
corrections of the energy of the ground and various excited states, we cal- 
culate the boundary critical exponents of single- and multicomponent Bethe 
ansatz soluble models. The boundary critical exponents are given in terms of 
the dressed charge matrix which has the same form as that of systems with 
periodic boundary conditions and is uniquely determined by the Bethe ansatz 
equations. A Luttinger liquid with open boundaries is the effective low-energy 
theory of these models. As applications of the theory, the Friedel oscillations 
due to the boundaries and the tunneling conductance through a barrier are 
also calculated. The tunneling conductance is determined by a nonuniversal 
boundary exponent which governs its power law dependence on temperature 
and frequency. 

PACS numbers: 05.70.jk, 64.60.Fr., 05.30.Fk 
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I. INTRODUCTION 



Phase transitions take place in a different way on surfaces and in tlie bulk of a sampleS. 
The exponents describing critical phenomena at surfaces differ from those of the bulk, and 
one may observe new phenomena due to anisotropy and the breaking of translational invari- 
ance caused by the boundary, like oscillations in correlation functions which, in the bulk, are 
monotonous, coordinate dependences, and in particular Friedel oscillations, in local quan- 
tities, etc. There are many approximate methods for calculating the critical properties of 
bulk and surface phenomena. 

Systems at a bulk critical point are not only scale invariant but also conformally invariant, 
a consequence of the combined rotational and scale invarianceHi. In two space dimensions or 
one space and one time dimension [thus including one-dimensional (ID) quantum systems], 
the constraints imposed by conformal invariance are much stronger than in higher dimensions 
because the conformal group is infinite-dimensional, and these constraints strongly simplify 
the calculation of correlation functions. The conformal field theory is parametrized by 
a unique constant - the conformal anomaly or the central charge c of the corresponding 
Virasoro algebraS. All critical exponents of theories with c < 1 are universal (independent 
of the interaction) and can be calculated exactlyi. However, for c > 1, the conformal 
dimensions continuously depend on the coupling of the fields and there is no general theory 
to deduce the critical exponents exactly. Conformal field theory only predicts that the 
central charge and the conformal dimensions can be derived from the finite size corrections 
to the energy and momentum spectrai'ii. For c > 1, it does not determine their actual 
numerical values. Still, this constitutes a powerful method to calculate the critical exponents 
of integrable ID quantum models because their energy and momentum spectra are known 
exactly, and this strategy has been applied successfully for bulk propertiesi0. 

An entirely parallel and equivalent development has taken place in the theory of correlated 
fermions in ID, using bosonization techniques and running under the name of "Luttinger 
liquid" The main focus here was on non- Fermi liquid properties of strongly correlated 
fermions, an exciting topic of current research. Fermi liquid theory fails in ID, and the 
Luttinger liquid provides the universal low-energy theory for gapless ID quantum systems. 
Its salient properties are: (i) absence of fermionic quasi-particle excitations, (ii) anoma- 
lous dimensions of various operators leading to non-universal power-law decay of correlation 
functions, and (iii) charge-spin separation. In terms of critical phenomena, the system is 
at a T = 0-quantum critical point, is conformally invariant, and the central charge of the 
associated Virasoro algebra(s) is unity. The exponents of the correlation functions (critical 
exponents) are related to each other by scaling relations and depend on a single renormal- 
ized coupling constant K per degree of freedom, playing the role of the Landau parameters 
familiar from Fermi liquid theory. For the Luttinger liquid, there are constitutive relations 
between three velocities characterizing the low-energy sector of the spectrum of the Hamil- 
tonian which determine the renormalized coupling constant and thus the critical exponents. 
For integrable models, the velocities and coupling constants have been determined from 
Bethe Ansatz (or other) solution^l. In some cases, conformal invariance has been used 
explicitly to determine the critical exponentsEi. For non-integrable models, they can be 
obtained reliably by exact diagonalization of the Hamiltonian. 

Recently, the problems of the one- dimensional systems with open boundaries have drawn 



2 



much attention. There are several issues of both experimental and theoretical relevance, (i) 
As in higher-dimensional systems, the boundary critical exponents are expected to be differ- 
ent from the bulk ones§. (ii) Experimentally studied systems are finite, and the progress in 
microfabrication of semiconductor structures has provided us with quantum wires so small 
that the boundary effects could become relevantEj. (iii) Here but also in the much bigger 
samples of quasi-lD metals, impurities can be relevant perturbations and effectively cut the 
systems into isolated strands of finite length. Specifically, a renormalization group analy- 
sis has shown that, for effectively repulsive interactions, the scattering potential due to an 
isolated impurity scales to infinitity in the low energy regime and thus the problem is equiv- 
alent to an open boundary proble An attempt to test experimentally the predictions 
of this theory has been published recently^. Such effects have also been invoked in the 
interpretation of electron spin diffusionlll and photoemissionllli§ experiments on quasi- ID 
organic conductors, (iv) Finally, some numerical methods such as the density-matrix renor- 
malization grou rely on the use of open boundaries and their results could, in principle, 
be affected by the chain ends. 



Both conformal field theory£3 and, very recently, the theory of Luttinger liquid*^, have 
been extended to systems with boundaries. What is missing to date, however, is an exact 
derivation of the boundary critical exponents from the Bethe ansatz solution of integrable 
quantum systems bridging the gap between microscopic (often lattice) models containing 
both high and low-energy physics, and the more effective theories for the low energy proper- 
ties. This gap has been bridged successfully for the periodic system It is the purpose 
of this paper to present such an exact derivation of boundary critical exponents. Moreover, 
in the course of the study, we shall see that the same strategy can be applied to determine 
the exponents of nonintegrable systems by exact numerical diagonalization provided they 
satisfy the basic assumption of conformal invariance. 

In this paper, we apply the method of the conformal field theory to the Bethe ansatz 
soluble models with open boundaries. The layout of the present paper is the following. In the 
following Section we briefly summarize some important results of boundary conformal field 
theory and of Luttinger liquid theory on bounded systems, in order to make the presentation 
self-contained and provide the basic tools. These are essentially the finite size corrections of 
the energy spectrum. In Section III we give a detailed calculation of the boundary critical 
exponents of two paradigmatic single-component Bethe ansatz soluble models (5-potential 
Bose gas and the antiferromagnetic Heisenberg chain). Section IV digresses to two important 
physical applications: the Friedel oscillation of the density distribution around the boundary 
(impurity) and the tunneling conductance through a barrier in a quantum wire. In Section V 
we generalize the result to the multicomponent case with the Hubbard model as an example. 
The summary in Section VI attempts to provide a broader perspective on our results. 

II. CONFORMAL FIELD THEORY AND LUTTINGER LIQUIDS IN BOUNDED 

SYSTEMS 

Systems with open boundaries 




^(x = 0) = ip{x = L) = 



(1) 
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(where is the wavefunction) are no longer space translational invariant but the time 
translational invariance is preserved. A two-point correlation function of a (primary) field 
at criticality then takes the general formS 

Gb{T, xi, X2) = (xiX2)"'^${[f + xl + xl]/xiX2} . (2) 

Eq. (0) applies to a ID quantum system, and r represents the imaginary time and Xi^2 the 
spatial coordinates; d is the conformal dimension of the primary field in the bulk; v is the 
Fermi velocity. When xi and X2 are near the boundary, and r ^ 00, Gb must behave as 

1 

Gfc(r,xi,X2) (3) 

2xb is the boundary critical exponent. Eq. (^ directly implies that limj,_>oo $(y) ~ y~^''- 
On the other hand, Cardy also showed that the ra-point correlation function in a half plane 
(with one open boundary at a; = 0) is identical to the 2n-point correlation function in the 
whole plane, provided only the ^-dependent part is taken into account in the latter. In this 
way, the two-point correlation function can also be represented as 

GbiZi,Z2,Zi,Z2) = [{Zi - Zi){z2 - Z2)l{Zi - Z2){Zx - Z2){zx - Z2)(^l " Z2)Y'^'^ Fb{y) , (4) 

where Fb{y) is an unknown scaling function. Here, we have switched to a notation in terms 
of complex variables Zj = vtj + ixj, Zj = vtj — ixj and y is given hy y = {zi — Z2){zi — 
Z2)l{z\ — z\){z2 — Z2)- For ?/ — > cxD, Fb{y) — > Direct comparison to Eq. (0) gives 

Xb = 4:d + a . (5) 

In the following text, we shall use Eq. (0) and Eq. (|^) alternatively. 

The conformal dimensions or critical exponents can be calculated from the finite size 
corrections of the energy spectra. To see this, consider the transformation 

C=—lnz, ( = —Inz , (6) 

TT 71 

applied to the upper half-plane a; > only. Such a conformal transformation maps the 
system from the semi-infinite plane onto a strip of width L with open boundary condition^ . 
From the general transformation properties of the correlation functions of a (primary) con- 
formal field (j){z, z) 

(0(.„.-O0(.2,^2)) = (||)^||)^(||)^(f|)^('^(Cl,Cl)0(C2,C2)) , (7) 

where A + A = c? and (j) and A — A = s give, respectively, the conformal dimension and spin 
of the field, one deduces for the correlation function on the strip 

JL sinh sinh ""^^^"^^^ 

(0(.„.-O0(z2,.2)). = Sinh ^ sinh ^ sinh ^II^M 

^ p sinh sinh 

sinh ^^^i^ sinh ^ ■ ^' 



There is also a spectral representation of the correlation function on the strip 



(0(2:1, ^i)0(z2,^2))l = Y.i^ix^,0)\n){n\<Pix2,0))exp[-ir^ - T2){EI - El)] , (9) 

n 

where E^ is the energy of the ground state and i?^^. the energies of the excited states; 
In) are the exact eigenstates of the Hamiltonian under consideration which form a complete 
set. Suppose E\ (the energy of the first excited state with the form factor (0|1) 7^ 0) takes 
the form (L — » 00), 

TTl) 1 

El-El = -x, + o{-). (10) 

From the asymptotic form of the correlation functions Eq. (^ and comparing Eqs. (^ and 
(^) we get Eq. (|^). Thus the finite size asymptotics of the low lying levels determines the 
boundary critical exponents. 

In general, the correlation function can oscillate, so that its asymptotics is not conformally 
invariant. In that case, however, one can decompose the field 0(^ z) into a sum of conformal 
fields 0„(z,z) which then determine the power-law asymptoticso. Because of the reflection 
symmetry of the open boundary systems, the field must have definite parity. We can 
expand odd and even parity fields as 

(j){z, z) = ^ (j)n{z, z) sminkpx) , ^(z, z) = ^ (pniz, z) cos{nkFx) , (11) 

n n 

respectively, n is an odd (even) integer for odd (even) parity fields. 

For systems with periodic boundary conditions, the Luttinger liquid phenomenology0lli 
provides a framework completely equivalent to conformal field theory but closer to the 
language of conventional solid state physics. It is based on the exactly solvable Luttinger 
model, and all physical properties can be described in terms of two parameters per degree 
of freedom [u = p,a for charge and spin), a renormalized sound velocity Vi, and an effective 
coupling constant Ki, which determines the decay of all correlation functions and thus the 
critical exponents. These parameters can be determined from the energies of the low-lying 
excited states of the Hamiltonianlll. 

This picture has been extended recently to systems with open boundarieJll. Due to 
the boundary conditions (|1|), the right- and left-moving fermions commonly used in the 
Luttinger model are not independent, and a single species moving, say, to the right [\l/+,s(x)] 
is sufficient, and it is periodic on a length 2L. The Fermi surface reduces to a single point 
+kp but the wavevectors k = mir/L > are quantized with twice the density of the periodic 
system. We then can rewrite the fermionic Hamiltonian 

Ho = -iVF J2 r «^,s(a;)5x^+,.(x) (12) 

in an equivalent form involving the bosonic density fluctuations (particle-hole excitations) 
p+^s{x) and "charge excitations" ANs corresponding to the addition of particles of spin s to 
the reference Fermi sea (i.e. Akps = ANsir/L) 
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Ho = 7rvpJ2 j dx : p+,,(x)p+,,(x) : +^5] (AiV,)^ . (13) 

The Fourier transform p^^s{p) of the density operators do not contain the p = 0-component 
which is represented exphcitly by ANg = Sfc(c+,s,fcC+,s,A; — (c+^s,fc'^+,s,A;)o) where (...)o denotes 
the (infinite) expectation value in the reference Fermi sea given by kp. Unhke the peri- 
odic case, "current excitations" describing the difference of right- and left-moving fermion 
numbers, cannot be defined in the bounded system. The Hamiltonian including forward 
scattering can then be diagonalized by a Bogoliubov tranformation as in the periodic case, 
defining the renormalized velocities Vi, of the bosonic charge and spin density fluctuations 
p+(p)[cr+(p)] = [p+^|(p) ± p+^|(p)]/v^, and coupling constants K^. The renormalized veloc- 
ity of the charge excitations ANp(^„-^ = AN-^ ± AA^^ is given by v^/ K^. The bosonization of 
this model is completed by an explicit representation of the Fermi operator \&_|_ s(x) in terms 
of the bosons p+ s(p)lli which allows to calculate all correlation functions of this model in 
terms of the Vy and Ky and thus defines its critical exponents. 

Vu and Ki, can now be found along the same lines as in the periodic systems^: (i) Vy 
can be computed directly from the spectrum of low-lying excitations; (ii) to get Kp, one 
calculates the compressibility k 

1 ^ 1 d'Eojn) ^ 2Kp _ 
K L dn^ TTVp 

The first equality gives the definition as the second derivative of the ground state energy 
with respect to particle density n = N/L and can be computed in the integrable model, 
and the second equality gives the Luttinger liquid expression which can be solved for Kp. 
Kcr is required to be unity by spin-rotation invariance, but can be calculated in the same 
way in more general cases from the susceptibility. Our determination of boundary critical 
exponents below can also be viewed as exploring this strategy. 



III. BETHE ANSATZ SOLUBLE MODELS 

We now compute the exact boundary critical exponents for Bethe ansatz soluble models 
with open boundaries. We first consider single- component models before turning to multi- 
component systems. In both cases, the critical exponents can be calculated explicitly. Two 
typical single- component models are the one dimensional ^-potential Bose gas model and 
the Heisenberg antiferromagnetic chain. Their Hamiltonians are 

Hbg = £ (<9^^^9^^ + c*^*'^^^ - /i^^^) dx, c> 0, h> 0, (15) 

Hxxz = eVj^ ■+! + c^J^jVi + cos 2r^a]cr]^, - ha]), (16) 
< 2?7 < TT, < /i < 4(1 - cos 2r]), 

where h is the chemical potential for the Bose gas model and the magnetic field for the 
Heisenberg chain. The anisotropy of the exchange integrals for the Heisenberg model is 
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Jz = cos(2?7); critical behavior is obtained only for easy-plane- type anisotropy \Jz\ < 1, and 
our definition of restricts us to this range, r] plays the role of a coupling constant, as is 
c for the Bose gas. a^'^'^ are the Pauli matrices. The open boundary conditions (|I]), for the 
Bose gas, translate into ^1/(0) = "^{L) = in terms of the boson operators "^{x), while they 
are contained in our representation (^) of the Heisenberg chain because the sites 1 and 
only couple to a single neighbor. 

These models are solved by means of the Bethe ansatzil. The A^-particle wave function 
is parametrized by numbers Xj which satisfy the equations 

2Lpo{X,) = 27r/, - 2y.(A,) - - Xi) + <I>(A, + A^)], (17) 

where po is the bare momentum and $ is the bare scattering phase: 

cosh(A + zr]) 

$^^(A) = -2arctan-, (18) 
smh(A — Ztr]) 

The numbers Ij are positive integers, and the parity effects known from periodic systems 
are absent in models with open boundaries. The bare energy of each particle is 

e^G(A) = X^-h, e^^^(A) = h-2 sinh^ 2r]/ cosh(A + ir]) cosh(A - ir]) . (19) 

The phase '^{Xj) = for the Bose gas and f{Xj) = Po{Xj) + $(2Aj) for the XXZ-chain. The 
eigenvalue of the Hamiltonian is equal to the sum of the bare energies of the particles 

N 

^L = Eeo(A,). (20) 

i=i 

Carefully checking the wave functions with the boundary conditions Eq. (|1|) we find that 
{±A} correspond to the same state. This is not surprising because of the reflection symmetry 
of the system, and corresponds to standing-wave-like solutions. The Bethe ansatz equation 
(p!7|) thus only allows solutions with Xj 7^ ±A/ for j ^ I so that all Xj > 0. The system 
therefore has only one Fermi point kp- This is very different from the case of the cyclic 
systems. A similar feature, however, occurs in Luttinger liquids in bounded systems (Section 
II). The appearence of a single Fermi point identifies the system as chiral, and ID quantum 
systems with boundaries therefore appear to be special cases of "chiral Luttinger liquids" , a 
notion that has appeared previously in the superficially quite unrelated area of edge states 
in the fractional quantum Hall effect^. 

With open boundaries, the total momentum is no longer a good quantum number. How- 
ever, the quantity 

N 

P = tJ:\I^ (21) 
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is conserved. We shall call it the "momentum" of the models. In the ground state, Ij takes 
consecutive integer values from 1 to (-^i=0 is not allowed). We can then define a dressed 
momentum 

p,(A,) = po(A,) + ^ Y.mX, - \i) + $(A, + A,)] + ^ = !^ . (22) 

The Fermi momentum kp = txN/L = 7m has the same value as in the periodic system. 

In the thermodynamic limit L —>■ oo, N oo keeping n finite, the ground state solves 
the following integral equations 

Pb{\) = ^ + r K,i\, ix)pb{ix)dfi , (23) 
vr JO 

rA 

eb(A) = eo(A) + / K^iX, p)ebip)dfx, (24) 
Jo 

Zn OA 

The cutoff parameter A is defined by the requirement e;,(A) = 0; Pb{X) is the density of A 
per unit length. Eq. (|2^) becomes 

PbiX) = poiX) + ^ /^[^(A -12) + ^(A + fi)]pbip)dfi . (25) 
2 Jo 

p'hiX) = Tiph{X) follows from a combination of (^) and (pSf), and ^^(A) = kp ensures that 
the Fermi surface is the same as in the finite system with the same electron density. 

We now compare equivalent quantities of the bounded and periodic systems and use the 
same symbols but without the subscript "6" for the latter. Directly comparing Eqs. (|23|)-(p5D 
with their periodic equivalentsi we find 

p,(A) = 2p(A) , 
e,(A) = e(A) , (26) 
Pb{X) =p{X) , 

a n = Uh (thus Af, = A). There is no restriction on sign(A) for the periodic systems but 
A > for the bounded ones. In the vicinity of the Fermi surface, the excitation energy is 

e{X)=v\pi,{X)-kF\, (27) 

where the Fermi velocity is given by 

e'(A) 4(A) 



p'{A) np,{Ay 



(28) 



It takes the same value of that of the system with the periodic boundary condition as it 
must be for the Fermi velocity can be determined by the leading term of the free energy 
which should not depend on the boundary conditions in the thermodynamic limit. 
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Unlike the cyclic systems, the systems with open boundaries have only two kinds of ele- 
mentary excitations: (i) Particle-hole (sound-like) excitations at the Fermi point kp. Their 
finite size corrections give the boundary critical exponent of the current-current correlation 
function for the Bose gas model and or the boundary critical exponent of the S'^-component 
of the spin correlation function for the Heisenberg spin chain, (ii) The change of the free 
energy induced by the variation of the particle number (termed "charge excitations" above). 
Its finite size correction gives the boundary critical exponent of the single particle correla- 
tion function for the Bose gas model and the critical exponent of the transverse spin-spin 
correlation function of the Heisenberg chain. These features are reproduced precisely in the 
Luttinger liquid theory of systems with open boundary conditions. 

We consider first the particle-hole excitations. To construct the lowest excitation state, 
we must put /at — Ja? + 1 in (|1^ keeping the other Ij unchanged. The change of the 
momentum is thus 



AP = ^, (29) 



and the excitation energy is 



AE,,L = 4 = 1. (30) 

To obtain the second kind of excitations (charge excitations), we define the dressed charge 
function as 

Z,(A) = 1 + r K{X, fi)Z,{fi)dfi. (31) 
Jo 

Obviously, Zb{X) = Z{X) for A > 0. The change of the free energy by AA^ additional 
particles isi 

A^,, . i|;„,„ , ^, _ -/.,„)] ^ (32) 
where /o is thefree energy density of the ground state. This gives 

^^''"^ ~ 2LZ2(A) ' - JZ^- ^^^^ 

For the general case, suppose a conformal field induce the momentum shift relative to 
the ground state as vrAJ/L and the change of the particle number AA^. The energy change 
relative to the ground state is then 

The boundary critical exponent takes the form 

2x,{AI, AN) = 2AI + (35) 
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where AI and AA^ are non-negative integers. The above relation was suggested by Alcaraz 
et a/.ii for integrable spin chains from numerical simulations. 

For example, we consider the Heisenberg chain in a zero magnetic field. In this casell, 
A = oo and 2Z^(oo) = tt/2t]. For ni, n2 near the boundary and t — >■ oo, the correlation 
functions take the following asymptotic forms 

«W<(0))~^, 

AJ = 1, AA^ = 0, r/ii =2x1 = 2; (36) 

KW<(0))~^, 
AJ = 0, AA^ = 1, ri^ = 2xi = ^ . (37) 

TT 

On the other hand, for t = 0, rii = 1 and n2 = n ^ 1, from Eq. @ we have 



K(0)<(0)) 



n 



0\\ = l + ^, for 2r^>- (38) 

^11 = 2 , for 2ri <- 
II , / - 2 



and 

K(0)a+(0)) 



0^ = ^ . (39) 



TT 



Notice that above we have used the bulk conformal dimensions d\\ = ^ for Arj > vr, = 1 
for ii] < TT and = ij/ii. When 2ri = n, the coupling is isotropic and r^y = rij_ = 2, 
^11 = = I as they should be for the spin SU{2) symmetry. The boundary critical 
exponents ?7||,_l measuring the decay of correlations with time, are two times larger than 
those of the bulk. On the other hand, there are new exponents 6\\^± for the decay of spatial 
correlations. In periodic systems, the two sets are identical as a consequence of conformal 
invariance. In bounded systems, the breaking of translational invariance along x by the 
boundary conditions while maintaining it along t, generates a new set of critical exponents. 
We now consider the finite size correction to the ground state energy 

N 

K,l = E^o(A[j7L]) . (40) 
i=i 

Using the Euler-Maclaurin formula we hav( 

El = Lj^e, (AM) dx + H- ^^U^. + ^^1^- + ^(Z) 

= LUn) + /, - - ^] + o(-). (41) 
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From Eqs. ([T9|) and (0) we know that efe(O) = 0. Thus we have 

Kl - Lfoin) ^h-^- (42) 

Here fb is the boundary energy which was extensively studiediiiiii'0. Eq. ( ^2] ) agrees with 
the predictions of conformal field theory!, and determines the central charge as c = 1. 



IV. FRIEDEL OSCILLATION AND TUNNELING CONDUCTANCE 

A. The Friedel oscillation 

Since the systems under consideration are obviously not translationally invariant, the 
density distribution is no longer homogeneous. Therfore, the ground state will exhibit 
Friedel oscillations. For a ID free fermion system of length L with open boundaries (p, the 
single-particle wave functions take the form 

vl/^(x) = y-sm-^, (43) 

with m positive integers. For the ground state, the density distribution in the box can be 
easily calculated as 

m=l m=l 

sin(2A;i;'x) r r i a a\ 

K,n ^ for X <t: L . (44) 

X 

For the interacting systems, we expect the density distribution to have a similar form 

{n{x))b^n ^— ^ —, 0<x<L, (45) 

where A and are two unknown constants and 7 is the exponent dominating the decay. 
As pointed out in Ref. 22, the n-point correlation functions of the open boundary systems 
are directly related to the 2n-point correlation functions of the periodic boundary systems. 
Thus 

{n{x))b = {n{z)n{z)) . (46) 

The right hand side of Eq. (|^) is to be understood in the sense that only the ^r-dependent 
part of the density-density correlation function (n(zi, Zi)n{z2, Z2)) is taken into account, and 
that we have to set Z2 = Zi. In addition, the oscillating term only originates from the current- 
current correlation {n{zi^Zi)n{z2.,Z2))j (in another languag , this object is designed as the 
2/ci?-charge density wave correlation function) which has been calculated in Ref. 9 as 

J^^i2kp{xi-X2) _j_ J^^i2kp{x2~xi) 

{n{z„ z,)n{z2, -Z2))j ^ - J^^^-^^^l^-—^^ . (47) 
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where S is a constant. Choosing the z-dependent part in Eq. p7|) and putting = -^i, we 
obtain Eq. (|45|) with 

7 = ^'(A) (48) 

This result quahtatively agrees with Ref. 15 where the Friedel oscillation for the Luttinger 
model with open boundaries was calculated, and determines the exact value of 7 from the 
Bethe ansatz equations through the dressed charge which can be calcualted easily. 



B. The tunneling conductance 

The boundary critical exponents are very important to study the tunneling effect in 
quantum wires. A strong barrier cuts the chain into two half-chains which, in first order, 
behave as two independent subsystems with an open boundary. For simplicity, we suppose 
them to be spinless and their Hamiltonian then reads 

^ = - E E { {cl,,Cr,j+i + H.c.) - Unr,jnr,,+i} - E E ^r,j, (49) 

r=± j=l r=±j=l 

where r = ± labels the two different half-chains; Clj {Crj) are the creation (annihilation) 
operators of the spinless fermions; rirj = CljCrj and h here denotes the chemical potential. 
The Hamiltonian is equivalent to an XX Z chain (|TB|) via a Jordan- Wigner transformation. 
We add a tunneling term to the Hamiltonian (|49|) 



T = -V 

The tunneling current is thus 



+ , y«l. (50) 



J=-iV 



- • (51) 



From linear response theory we know that the tunneling conductance up to order V'^ is given 
by 

G{uj) = % J dte'^' J dt'e{t - t') J dt"e{t - t"){[J{t), J{t")]) . (52) 

Since the average (■ ■ ■) is taken at T = 0, the current correlation function may be separated 
into 



m), Jit")]) ~ n {CrAt)cun) -it- fr^^- , (53) 



where ri±_ is the boundary critical exponent of the single particle correlation function. Sub- 
stituting the above relation into (E2|) we readily obtain 



Giu) r^J, e = 2iv^ - 1) = - 2 . (54) 
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At finite but very low temperatures T ~ 0, the conductance behaves as 

G(T) = GoT', (55) 

where Go is a constant. 

At U = 0, H describes free fermions with a barrier. In this case, r]± = 1 and G{T) is 
independent of temperature and finite. The system is marginaL For f/ > 0, 77^ > 1 and 
the conductance tends to zero as T — 0. The fermion-barrier scattering is relevant and the 
"Coulomb blockade" behavior arises - a result consistent with the observations of Kane and 
FisherEl. For U < 0, r]± < 1 and the tunneling conductance diverges as T — * 0. This is 
a consequence of the divergent superconducting fiuctuations found in that situation. Our 
calculation may be generalized to general Luttinger liquids and the critical exponent 6 is 
then given by 

9 = Y.K;'-2, (56) 

V 

where are the stiffness constants of the Luttinger liquid under consideration. 



V. MULTICOMPONENT INTEGRABLE MODELS 

Recently much attention has been focused on the open boundary problem for inte- 
grable models with mult i- component fieldsilli'i. Typical models are the one-dimensional 
5-potential Fermi gas modell^, the Hubbard chainE^I and the supersymmetric t — J model 
with open boundarieslil'ii. The above discussion can also be generalized to these models. In 
these cases, the refiection Bethe ansatz equations take the general form 

M Np 

2LK(AJ) = 27r/; - E E't-^^/slA- - Af) + ^./^(AJ + Af)], (57) 

/3=1 1=1 

where $0,^3 (A° — Af) are the bare scattering phases and odd functions of their arguments, 
M is the number of the components, and the prime after the sums means that when a = j3, 
j ^ l. The eigenvalue of the Hamiltonian is 

M Na 

^=EEeo(A,")- (58) 

a=lj=l 

Also, A" > is supposed. In complete analogy to Eqs. (|23|)-(|25|) and (|3TD, dressed quantities 
are defined as 



pUXn = ^-^^^ + E / ' ^"MA", A^)pf (A^)rfA^ , (59) 

TT Jo 

1 r^0 r la 

K(A") = Poin + 7;T. <f ./^(A" - A'^) + <^'a,(A" + A'^) p',{X')dX^ , (60) 

e°(A") = e^(A") + / ^./.(A", A^ef (A^)rfA^ , (61) 
,3=1 -^0 

ZipiX^) = 5a/5 + E / Zi^{\^)K,p{X\ A^)c?A^ , (62) 

1 



7=1 
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where /^^/^(A", A^) = ^[$„/3'(A" - A^) + $c./3'(A° + A^)] is an even function. The relations 

= vrp^(A") = 2vrp.(A") , (63) 
e^(A") = e.(A°) , (64) 
ZU>^') = Z^,{\^) , (65) 

for A" > 0, compare bounded and periodic systems at equal density and generalize (p6D. 
The finite size correction to the energies of the excited states is then 

E,,L - El, = ^^^^a |-[(Z-iAN)j2 + A/.j + o(-), (66) 

where Zap = ZapiAfs), AN = {ANi, ANm} are M-dimensional vectors with integer 
components. The number ANa gives the change of Na, the number of pseudoparticles of 
type a (pseudoparticles refers to the particle-like excitations in the Bethe ansatz and not 
necessarily to physical particles), in the excited state with respect to the ground state (i.e. 
the charge excitations). The non-negative integers Ala describe pseudoparticle-pseudohole 
excitations (more precisely, a change of X^j^iPb (A") in units of ^) in the vicinity kp (Fermi 
momenta kp of the pseudoparticles are defined as kp = p^{Aa) = una). The Fermi velocity 



is Va = ^" it \ as usual 



ea'(Aa) 

The finite size correction of the ground state energy is given by 

EIl = - E / "K'(A")enA")rfA" + A - ^ E ^« + ^(t)- (67) 



The Fermi velocities Va are arbitrary in principle and quantitatively depend on details of 
the interactions in practice. As a consequence, the system is described by a sum of M 
conformal algebras, each with a central charge 1. Their contributions to the boundary 
critical exponents 

= [(Z-i AN),]2 + 2AIa, (68) 
are additive, and the total boundary critical exponent is thus 

M M 

2xb = 2j2xb =2Y,^Ia + (Z-^AN)^(Z-^AN). (69) 

The same structure is found in periodic systems. 

As an example, we give some leading boundary critical exponents of the Hubbard chain 
with open boundaries. The bulk critical exponents of this model were determined by Frahm 
and Korepin0. The Hamiltonian reads 

N-l N N ^ TV 

^ = - E E c'/aC^+l<x + 4f/ E ^n^n - E E ^^'^ - 2 E("'iT - ^n)^ (7o) 

i=l (T=± i=l 1=1 (T=± i=l 

where Cio- (C|^) is the electron annihilation (creation) operator; /i denotes the chemical 
potential and h is the external magnetic field. The wave functions are parametrized by two 
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sets of parameters {k} and {A}, the rapidities of the charges and spins, respectively. The 
following set of integral equations determine their bare (e^^) and dressed (e^'"*) energies and 
distribution functions {pl''^) 

h h 
eo{k) = -2 cos k + fi--, eo(A)=2' 

e^(fc) =eS(fc)+ r^Ki(sinA;,A)e^(A)rfA, 
Jo 

e^(A) = e^(A)+ / ' cos A:/s:i(A, sin A;)e;;(A;)dA; - / V2(A, /i)e^(/i)d/i, (71) 
Jo Jo 

1 rAs 

PbW = - + COS / i^'i(sin k, \)pUX)dX, 

TT JO 

pt{X)= f^^ K^{X,smk)pl{k)dk- j^' K2{X,p)pl{p)dp. (72) 
Jo Jo 



The dressed charge (0), in the multicomponent problem, takes a matrix structure, with 
elements 

ZUk) = 1 + Zl{X)Ki{X,sink)dX, 
Jo 

Z^^XX)= I \oskZ^^^{k)Ki{smk,X)dk- [ ^ Zl{p)K2{p,X)dp, 
Jo Jo 

Zlik) = Zl{X)K,{X,smk)dX, (73) 
Jo 

Zl{X) = l+ [ \oskZ,cik)Ki{smk,X)dk- f ^ Zl{p)K2{p, X)dp, 
Jo Jo 

with the kernel 

1 TlU TlU 

= ^\nUy + (A - py + {nUr + (A + p^^^ " = ^^'^ 

As in the single-component case, the dressed charge matrices for open and periodic boundary 
conditions are identical, = Zq,^^. The Fermi velocities are given by 

and Ac,s are defined by e^''*(Ac,s) = 0. 

Below we list some correlation functions and the ANc^s, AIc^s which must be used in 
Eqs. ([68| ) and (|69D in order to determine the leading critical exponents: 

(1) The field correlator 

G^^ix^,X2,t) = (a,TWCl.T(0))' (76) 
AN, = 1, ANs = 0, A/, = A/, = 0. 

(2) The density-density correlator 



15 



Gnnixi,X2,t) = {n^^{t)n^^{0)) , (77) 

A/, = 1, A. 

(3) The spin-spin correlators 



AN^ = ANs = 0, A/, = 1, AIs = or AI^ = 0, AI^ = 1. 



GU^i,X2,t) = {S'ix,,t)S'ix2,0)) , (78) 
S"{x,t) = - [n^^ - n^J, 
ANc = ANs = 0, Ale = 1, AIs = or AI^ = 0, A/, = 1. 

GUx,,X2,t) = {S-{x^,t)S+{x2,0)) , (79) 

5+(x,t) = cl^(t)ai(t), 

AA/'c = 0, ANs = 1, A/c = AIs = 0. 

(4) The triplet pair correlator 

G^^\xuX2,t) = (a,+iT(t)a,T(i)C^l.T(0)Ci,+iT(0)) , (80) 
AN^ = 2, ANs = 0, A/e = A/, = 0. 

(5) The singlet pair correlator 

Gf\xuX2,t) = {C,At)C.At)CU0)Cl^m , (81) 
AA^, = 2, ANs = 1, A/e = A/, = 0. 

Precise values for the critical exponents then follow immediately, via (^) and (|69D, once the 
dressed charge matrix (^) is calculated. This is a matter of routine, and due to the equality 
of this matrix for open and periodic systems (cf. above), the published results for the periodic 
Hubbard mo del0 can be used directly to evaluate the boundary critical expoenents. 



VI. SUMMARY 



We have derived explicitly the boundary critical exponents of both single- and multi- 
component Bethe ansatz soluble models of interacting bosons and fermions. Our results 
imply that the descendant fields (particle-hole excitations) contribute the same (integer) 
amount to the boundary and the bulk critical exponents. However, the contribution from 
charge excitations (additional particles) to the boundary critical exponents is twice as big as 
to the bulk exponents. The current excitations are completely depressed for open boundaries 
and thus contribute nothing to the boundary critical exponents. Apparently, this statement 
is valid much beyond the the Bethe ansatz soluble models and applies in general to Luttinger 
liquids with open boundariesEl. The critical exponents are determined by the dressed charge 
matrix which we have shown to be independent of the boundary conditions. Moreover, our 
method of calculation relies only on the determination of energies which can be performed 
accurately by numerical methods in models which cannot be solved by Bethe ansatz. There- 
fore, one can determine, at least numerically, the boundary critical exponents for all ID 
quantum systems, provided they are conformally invariant, by the method described in this 
paper. 
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